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. Recent experiments on quantum behavior in microfabricated solid-state systems suggest tantaliz- 

^ ' ing connections to quantum optics. Several of these experiments address the prototypical problem 

04 I of cavity quantum electrodynamics: a two-level system coupled to a quantum harmonic oscillator. 

Such devices may allow the exploration of parameter regimes outside the near-resonance and weak- 
coupling assumptions of the ubiquitous rotating- wave approximation (RWA), necessitating other 
theoretical approaches. One such approach is an adiabatic approximation in the limit that the oscil- 
lator frequency is much larger than the characteristic frequency of the two-level system. A derivation 
of the approximation is presented, together with a discussion of its applicability in a system con- 
sisting of a Cooper-pair box coupled to a nanomechanical resonator. Within this approximation 
, the time evolution of the two-level-system occupation probability is calculated using both thermal- 

and coherent-state initial conditions for the oscillator, focusing particularly on collapse and revival 
phenomena. For thermal-state initial conditions parameter regimes are found in which collapse and 
revival regions may be clearly distinguished, unlike the erratic evolution of the thermal-state RWA 
, model. Coherent-state initial conditions lead to complex behavior, which exhibits sensitive depen- 

dence on the coupling strength and the initial amplitude of the oscillator state. One feature of the 
regime considered here is that closed-form evaluation of the time evolution may be carried out in the 
, weak-coupling limit, which provides insight into the differences between the thermal- and coherent- 

51^ ' state models. Finally, potential experimental observations in solid-state systems, particularly the 

O ' Cooper-pair box-nanomechanical resonator system, are discussed and found to be promising. 

PACS numbers: 42.50.Md, 42.50.Hz, 85.25.Cp, 85.85.-|-j 
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^ . I. INTRODUCTION 

o\ 

' One of the simplest fully quantum-mechanical systems consists of a harmonic oscillator coupled to a two-level (spin- 
" like) system. Introduced in 1963 by Jaynes and Cummingsji this model has yet to be completely solved despite its 
' apparent simplicity. The approach taken by the original authors, known as the rotating- wave approximation (RWA), 
\ relies upon the assumptions of near-resonance and weak coupling between the two systems. The RWA is widely used 
. because it is readily solvable and describes quite accurately the standard physical realization of such a system: an 
' atom coupled to a field mode of an electromagnetic cavity. In this experimental situation, the coupling strength 
, between the atom and the field is largely determined by the intrinsic dipole moment of the atom; for all experiments 
I ' to date, the coupling strength A is very small compared to the atomic transition frequency uj {X/uj ^ 10~^-10^®).^j^ 
' ^ The near-resonance condition is necessary to ensure validity of the two-level description of the atom. Thus the RWA 

is a natural, and excellent, approximation in such a system. 
^ Quantum-limited solid state devices offer an alternative to the traditional atom-cavity implementation of the spin- 

oscillator system. Recent experiments have shown clear spectroscopic evidence that a Cooper-pair box (CPB), or 
Josephson charge qubit, coupled to a superconducting transmission line behaves much like an atom in a cavity. The 
dipole coupling between the two systems is X/uj « 10""^, 3 — 4 orders of magnitude larger than that achieved in 
^ ' atomic systemsi^ Capacitive or inductive couplings offer the possibility of still larger coupling strengths than those 
5^ ! possible with dipole coupling, even at large detunings between the fundamental frequencies of the oscillator and the 
two-level system. Some results from a flux-based, inductively-coupled system give preliminary evidence for coupled 
quantum behavior and entanglement between the two-level system and the oscillator.^ Another intriguing possibility 
involves capacitively coupling a CPB^SiS or a Josephson phase qubit-'^-'i to a nanomechanical resonator (NR). All of 
these systems are capable of accessing coupling strengths and detunings outside the regime in which the RWA is valid, 
requiring different theoretical approaches to the problem. For example, Ref. ^3 details a perturbative treatment which 
is valid for moderately strong coupling {X/uj < 0.3) at zero detuning. 

In this paper we discuss an approximation which is able to treat strong coupling and large detuning. It is valid when 
the splitting frequency of the two-level system is much smaller than the frequency of the oscillator and holds well even 
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for coupling strengths up to or larger than the oscillator frequency. The approximation is used to examine the time 
evolution of the two-level system when the harmonic oscillator begins in a thermal state or a coherent state. Several 
effects of the coupling to the oscillator are distinguished, including enhanced apparent decoherence rates, frequency 
modification, and collapse and revivals of Rabi oscillations. We focus particularly on collapses and revivals in this 
model. 

In the next section we introduce the form of the Hamiltonian to be considered and derive the adiabatic approxima- 
tion. In Sec, mil we evaluate the approximate time evolution of the two-level system coupled to, respectively, a Fock 
state, a thermal state, and a coherent state of the oscillator, and classify the behavior in several parameter regimes. 
Section llVl contains a further approximation which allows evaluation in closed form of the infinite sums encountered in 
the thermal- and coherent-state models. The prospects for experimental observation of our predictions are analyzed 
in Sec.0 and Sec. IVII concludes the paper. The Appendix contains a detailed derivation of the effective Hamiltonian 
for a CPB-NR system. 



II. ADIABATIC APPROXIMATION IN THE DISPLACED OSCILLATOR BASIS 

The Hamiltonian which forms the basis for the calculations and discussion in this paper is 

H = ^hflcTx + hXazia^ + a) + huQa^a. (1) 

The Hamiltonian typically considered in cavity quantum electrodynamicsi^ii (CQED) differs from Eq. |^ only by a 
rotation about ay. Our notation is chosen based on the particular solid-state system which we have in mind, consisting 
of a Cooper-pair box coupled to a nanomechanical resonator. When the CPB is biased away from its degeneracy point 
and driven by a classical field resonant with the qubit transition frequency, an effective Hamiltonian for the coupled 
system may be found which has exactly the form of Eq. |^ (see Sec.lVland the Appendix) In this case the two- level 
system undergoes Rabi oscillations due to the (classical) external driving field which are then altered by the coupling 
to the quantum oscillator, unlike the usual picture in which the two-level system is driven by the quantum oscillator 
itself. Although the Hamiltonian and therefore the results we obtain are not restricted to the CPB-NR system, for 
concreteness we will use the language of that system throughout most of the paper. 

As no analytic solution to Eq. is yet known, some approximation is required. The most common approach 
is to assume that the two-level system and the oscillator are close to resonance, (fl — loq) <^ ^^,wo, and that the 
coupling between them is weak, A ^ ^,ujo- Then terms which do not conserve energy may be discarded (the RWA), 
resulting in what is usually termed the Jaynes-Cummings model (JCM)?^*^^ As mentioned in the Introduction, this 
is an appropriate approximation for atom-cavity experiments; however, unlike atom-cavity systems, the solid-state 
system we are considering has the potential for strong coupling at large detunings. In this paper we will treat the 
case in which the two- level-system splitting frequency is much smaller than the oscillator frequency, f2 <C wq, and the 
coupling strength is allowed to be large, on the order of or greater than the oscillator frequency. The RWA is not 
appropriate in this limit. 

An excellent way to treat this regime is via a type of adiabatic approximation. This approximation has been 
derived previously by several authors using different methods. Graham and Hohnerbach refer to this regime as 
the "quasi-degenerate limit" and give the same lowest-order expressions as we deriveii2iiii4 Schwcber utilized the 
Bargmann Hilbert-space representationi^ and Crisp solved recurrence relations;^® both of these authors found higher- 
order corrections beyond what we present. We take yet a different approach. First, by neglecting the self-energy of 
the two-level system, we derive the basis in which the rest of the calculation will be performed, called the "displaced 
oscillator" basis. In this basis the Hamiltonian may be truncated to a block diagonal form and the blocks solved 
individually. Essentially, the two- level- system self-energy is introduced only as needed to lift the degeneracy within 
individual subspaces in the displaced oscillator basis. 

To begin with we consider eigenstates of the form = \i) (8) \4>i) where \i = +, — ) denotes the eigenstates of 

az and are the corresponding oscillator eigenstates, found from the last two terms of Eq. ^ with az set to its 
eigenvalue of ±1 as appropriate: 

[±hX{a^ + a) + hujoa^ a]\(j)±) = E\(P±) . (2) 

Completing the square gives 




Taking X/ujq to be real, the operator on the left-hand side may be rewritten as 

^D{TX/uJo)a^aDHTX/uJo), (4) 
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FIG. 1: Graphical representation of the displaced oscillator basis. The equilibrium position of the harmonic oscillator well is 
shifted by an amount proportional to the coupling constant A, with the direction of the shift determined by the state of the 
two-level system. Each well retains its harmonic character, with the usual eigenstates. Eigenstates with the same value of A'' 
are degenerate in energy. 



where D{iy) = exp[i^(d^ — a)] is a displacement operatorjii The operator in Eq. 21may be interpreted as the number 
operator for a harmonic oscillator of frequency loq, effective mass m, and equilibrium position =p(2A/wo)Aa;zp where 
Aa;zp = \fhJ(^2MiL0Q) . The eigenstates of this operator are displaced Fock (number) states, 

|c/)±) = e=F(^/'^°)(^^-^)|A^) = |iV±) A^ = 0,l,2,... (5) 

with energies given by 

En ^hw^(N (6) 

These approximate eigenstates and energies constitute the displaced oscillator basis, which will be used throughout 
subsequent calculations. An illustration of the displaced- well potentials and the corresponding energy levels is given 
in Fig.n 

This basis has a simple interpretation in the CPB-NR system. Due to the capacitive coupling, one charge state of 
the CPB will attract the NR and shift its equilibrium position closer to the CPB, while the oppposite charge state will 
repel the NR. As long as the stretching of the beam due to the attraction is small, the NR will remain approximately 
linear and may still be described as a harmonic oscillator. However, the position of the harmonic potential well will 
have shifted. This is precisely the situation described mathematically by the displaced oscillator basis, which allows 
the use of the ordinary harmonic oscillator formalism within each displaced potential well. The correspondence of 
physical position displacement to the formal displacement operator makes this a natural basis to work in. 

However, there are some peculiarities associated with using the displaced oscillator states as a basis. Although the 
states |-|-, A^+) (| — , A^_)) form an orthonormal basis with (M+|iV+) — Smn {{M-\N-) = 5mn), the states and 
|A^_) are not mutually orthogonal. This may be understood easily in position space, where the displacement operator 
corresponds to a displacement in x. Consider, for example, the harmonic oscillator ground state |0), which has a 
Gaussian form: this state displaced by a finite amount is never completely orthogonal to the same state displaced by 
the opposite amount, although the displaced states may become very nearly orthogonal if the displacement is large 
enough. The overlap between Fock states displaced in different directions is given by 




/M\N\^r • '■v-2AM)*^-^yM7M!<-^[(2A/a;o)^] M>N 



where Ll is an associated Laguerre polynomial. Note that (Af_|iV+) = (—1)^ *^(A^_|Af+). Together with the fact 
that {M_\N+) is real, this imphes that {M+\N_) = (-l)^^~^(M_|iV+), which is a useful identity The lack of 
orthogonality between different displacements leads to the unusual results in the two-level-system dynamics which 
will be found later. 

The next step is to reintroduce the two-level-system self-energy, 17, which is assumed to be small. Reformulating 
the problem in terms of a matrix written in the displaced oscillator basis allows the approximation to be carried out 
in a natural way. Using the overlap functions calculated above, the matrix representing the Hamiltonian may be 
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FIG. 2: Graphical representation of the adiabatic approximation. The two displaced oscillator wells, illustrated in Fig. are 
now allowed to interact. Levels with the same value of A'' become mixed; the resultant energy splitting is proportional to the 
overlap of the wavefunctions. 
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In this approximation, the self-energy of the two-level system has been employed only where it is needed to lift the 
degeneracy of the displaced-oscillator basis states [Eqs. lO and lO; also see Fig.^. The natural interpretation of this 
form of the Hamiltonian is that the oscillator will be restricted to remain in the N± subspace if it has been initialized 
in |A'^+) or \N-). In other words, mixing occurs only between levels in opposite wells which have the same value of N 
and thus the same energy. In this picture the condition f2 <C may be seen as a statement about statics: the spacing 
of the oscillator energy levels is very large compared to the spacing of the two-level system, so a transition in the 
two-level system can never excite the oscillator. An alternative argument, found in Sec. III(C) of Ref.IIl relies on the 
separation of characteristic times in the two systems: given the assumption Q. <^ ojq, the oscillator responds almost 
instantaneously to changes in (o-^), so that the oscillator dynamics is "slaved" to the two-level-system dynamics. 
Due to its simple block diagonal form, Eq. @ may be solved easily. The solutions in the A^th block are given by 



\^±,n) = ^{\ + ,N+)±\-,N^)1 
E±^N ^±]p{N^\N+)+EN. 



(10a) 
(10b) 



These energies and eigenstates constitute the adiabatic approximation to lowest order in fi/wo- A graphical represen- 
tation is shown in Fig. |2| 

Figure l^l plots the lowest-lying energy levels given in Eq. (jlObp as a function of the coupling strength A/wq for 
different values of the ratio ^/loq. For comparison purposes, the results of a numerical diagonalization of the full 
Hamiltonian [Eq. are also shown. 
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FIG. 3: (Color online) Energy levels given by analytic approximation methods (solid lines) and by numerical solution of the full 
Hamiltonian (dashed lines), (a) Adiabatic approximation and numerical solution with Q,/loq — 1/3. (b) Adiabatic approximation 
and numerical solution with Q/ujo ~ (c) Rotating wave approximation and numerical solution with Q/uio = 1. (d) Adiabatic 
approximation and numerical solution with Q/ujo — 3. 

Consider first the case H/ojo = 1/3, shown in Fig.lSfa). This is within the regime in which the approximate solution 
should be vaUd. The "ripple" structure imposed by the Laguerre polynomials on the smooth variation of E 
with X/uJo is immediately evident. This structure was noted in Ref. ^3 3.nd interpreted as an interference between 
states displaced in opposite directions. In the limit X/loq —^ oo the distance between the wells becomes infinite, and 
the overlap (7V_|iV+) — > 0. With no mixing between the wells the spectrum becomes that of two identical harmonic 
oscillators and the energy levels become pairwise degenerate. Agreement with the numerical solution is excellent. 

Another noteworthy feature of the structure of the energy levels is the multiple crossings which appear between 
pairs of levels. These are true crossings, not narrow avoided crossings, allowed by conservation of the parity operator 
P — exp[i7r(a^a + 5 + ^o-^)] in Eq. (Qii^ The approximate eigenstates HlOaf) are eigenstates of P; pairs of levels with 

different eigenvalues of P are allowed to cross. In other words, at the crossing points the quantum states in the two 
displaced wells destructively interfere with one another, destroying the tunneling process which mixes the states and 
provides the energy splitting. 

Figure |2tb) shows the resonance case, il/ujo = 1. Plots for this case also appear in Ref. H^i- Interestingly, the 
agreement between the approximate solution and the numerical solution is still quite good, especially at larger coupling 
strengths. Compare this to the RWA energy levelsfi = NhujQ ± n/2 ± HXVN, shown in Fig.EIc). The RWA 

gives the correct limiting behavior as X/loq — > 0, but diverges from the numerical solution starting around the point 
where the paired levels first cross. This comparison illustrates the dependence of the RWA upon the assumption of 
weak coupling, even when the oscillator and two-level system are exactly resonant. 

Finally, Fig. |2[d) shows the case fl/ujo — 3, which is out of the regime in which the adiabatic approximation is 
expected to hold. Plots of this case are given in Ref. 0; however, the approximate and numerical solutions are shown 
in different figures, making it difficult to compare the two. The plot shown here demonstrates that the adiabatic 
approximation breaks down in this regime except in the broadest qualitative sense. When fl/ujo > 1 spurious level 
crossings appear at small values of X/ujq, and a substantial phase difference develops between the ripples in the 
numerical and approximate solutions. Thus it may be seen that, as expected, the adiabatic approximation is not a 
reasonable treatment for Q/loq > 1. 

The displaced oscillator basis provides a physically intuitive picture for the derivation of an adiabatic approximation 
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in the regime rj/wo < 1. Comparing the adiabatic approximation with a numerical solution of the full Hamiltonian 
shows that the approximation works quite well in the regime for which it was derived. Although this approximation 
has been derived previously, no authors seem to have explored its consequences for experimental observables. Much 
of the remainder of this paper will be devoted to a study of those consequences. 

III. DYNAMICAL BEHAVIOR 

For applications in real systems, the dynamical behavior of the two-level system is of particular interest. In this 
section we discuss the time dependence of the two- level-system observable ct^, which corresponds to charge in a 
Cooper-pair box. Three commonly used states are considered as initial conditions for the harmonic oscillator: the 
Fock state, the thermal state, and the coherent state, each of which may be applicable in different situations. As 
might be expected, the behavior of the two-level system changes dramatically depending on the initial oscillator 
state. We have verified the behavior obtained in the adiabatic approximation against a numerical solution of the full 
Hamiltonian. Provided the ratio uj^^/^l is made large enough, the agreement is excellent over the short time scales 
which are likely to be experimentally accessible; the value ojq/^ = 10 gives quantitative agreement over several dozen 
periods of oscillation. Time evolution in the adiabatic approximation shows a rich variety of behavior which we 
demonstrate and classify. 

Throughout this section the quantity to be examined is the probability of obtaining the state |— ) as a function of 
time, P(— , t). The initial state of the two- level system is taken to be |— ), and the initial state of the oscillator is given 
in the displaced basis corresponding to the state |— ). This situation might be obtained, for example, by turning on 
the bias voltage Vg between the NR and the CPB, tuning the CPB gate bias voltage Vb so that the net bias is away 
from the degeneracy point, and allowing the CPB to relax to its ground state. Preparation of the oscillator state 
would vary depending on the type of state desired (Fock, thermal, or coherent); the discussion of each state includes 
some indication of how that state might be prepared. At time t = Q the microwave field is switched on and the system 
begins to evolve in time. Other initialization schemes might be imagined, but starting the oscillator in the displaced 
basis simplifies the mathematics. 

An important point to keep in mind when comparing the results presented here with results from the JCM is 
that they are measured in different bases relative to their respective Hamiltonians. In our notation, the state |— ) is 
an eigenstate of the two-level-system operator az which is associated with the interaction term of the Hamiltonian. 
In the JCM the initial state and the measured state of the atom are typically chosen to be eigenstates of the bare 
atomic Hamiltonian rather than the interaction Hamiltonian. These are natural choices for the corresponding physical 
systems, but it should be kept in mind that a degree of caution must be used in comparing the two models. 

The simplest dynamical behavior is obtained when the oscillator begins in a displaced Fock state, such that |'!/'(0)) = 
|— ) (X) \NJ). The time evolution of the probability to find the two-level system in the initial state |— ) is given by 

p{-,t) = \{-,N^\m)? 

= cos2(17'Wt/2), ^ ' 

where 

= n{N^\N+) = f^e-2-^'/'^«ijv[(2A/cjo)^]. (12) 

The CPB undergoes Rabi oscillations (due to the external driving field), with a frequency which is modified by the 
strength of the coupling to the resonator; following Ref. ^| we refer to this as "adiabatic renormalization."— Fig. ^ 
shows a plot of the renormalized frequencies H12|) versus the coupling strength A/wq for small values of N . Unlike 
the Rabi frequencies obtained in the JCM, the frequencies found here are not monotonic functions of the coupling 
strength A or of for > 0. 

One unexpected feature of these frequencies is the zeros which occur as the coupling strength increases. Comparison 
with a numerical solution of the full Hamiltonian shows that the adiabatic approximation breaks down to some extent 
around the "critical points" in the coupling strength at which the renormalized Rabi frequency goes to zero. This 
point is discussed in Ref. J£j and a higher-order formula is derived which is valid even near the critical points. Some 
caution must therefore be used in making predictions for Fock state initial conditions from our lowest-order formula. 
However, Fock states are highly non-classical states and although some methods for preparing such states in NRs have 
been proposed f^ii^i the experiments appear difficult. This paper is primarily concerned with time evolution from more 
realistic initial states which involve some distribution of number states. The distribution reduces the contribution 
from any given number state, and our numerical studies indicate that the approximation works well even when states 
are included which have critical points near a given coupling strength. 



7 




The next initial condition for the oscillator which we consider is a thermally-occupied state, sometimes also referred 
to as a chaotic state, in the displaced basis. This type of state is expected for an oscillator in thermal equilibrium 
with its environment, such as the NRs studied in the recent experiments of LaHaye et al.^*^ The thermal state is a 
fully mixed state which must be described by a density matrix rather than a state vector. The two-level system is 
again taken to be initialized in the state |— ), so that the initial density matrix for the coupled system is given by 



Pth(O) = !-)(-! ®5]pth(^)|A^-)(A^-l, (13) 

where 



N 



^''^^''^ ^ (l + (Ar))(l + l/(iV))iv (14) 

and (N) = {gfi-^o/kBT _ ^^-i jg ^j^g average number of quanta in the oscillator at the temperature T. 

Assuming that the system is weakly coupled to the thermal environment so that the influence of the environment 
is restricted to determining {N) (i.e., decoherence is not included), the time evolution of the system is given by 

Pth(i) = e-*^*/Vth(0)e^^*/'' 

= ^pth(A^)[cos2(f7'(^)t/2) I-, 7V_) (-, 7V_| + sin2(f7'(^)V2) |+, N+) (+, N+\ 

JV 

+ ^ sin(f^'Wt/2) cos(l]'Wt/2) (|-, iV_)(+, 7V+| - | + , iV+)(-, 7V_|)]. 
The reduced density matrix for the CPB is obtained by tracing over the oscillator states: 

PtT^t)=E^''h(iV)[cos2(17'Wt/2)|->(-|+sin2(J7'Wt/2) |+)(+| 

JV (15) 
+ isin(r!'Wi/2) cos(r!'Wt/2) (A^-|iV+) (|-)(+| - |+)(-|)]. 

From the reduced density matrix the probability of obtaining the state | — ) is found to be 

Pth(-,t) = (-|p?r(t)|-) 

JV 

If no further approximations are made, the sum in Eq. requires numerical evaluation. The parameter space is 
complicated, but at least three qualitatively distinct regimes of behavior may be found, characterized by the oscillator 
temperature and the coupling strength. In the very low temperature regime, {N) k, 0.01, the behavior consists of 
ordinary Rabi oscillations with a frequency renormalized by the coupling to the ground state of the oscillator. This 
rcnormalization becomes significant for relatively large coupling strengths, X/loq « 0.1-1. 

However, due to the large width in N of the thermal distribution, the different frequencies fi''^^ involved in the time 
series ()16|l tend to interfere with each other, producing a decay in the amplitude of the oscillations. This effect, known 
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FIG. 5: (Color online) Behavior of the two-level-system occupation probability in the high-temperature ((A^) = 100), weak- 
coupling regime for short times (main figure) and long times (inset). Throughout the paper all time-dependent plots have 
Q./LJO = 1/10. 




FIG. 6: (Color online) Behavior of the first revival region as a function of temperature, with X/uJo = 0.1. 



as a "collapse" in quantum optics, ^^'^^ dominates the short-time behavior at higher temperatures as long as X/coq 
is not too large. At longer times, the discrete nature of the spectrum allows a partial rephasing of the oscillations, 
resulting in "revivals" of the oscillation amplitude. These phenomena are illustrated in Fig. |S1 Several features are 
worth noting. First is the nonzero amplitude in the collapse region; this residual amplitude decreases with increasing 
temperature. Second, the collapse time and the times at which the revivals occur depend on the coupling strength: 
both the collapse time and the revival time shorten as the coupling strength increases. Finally, the amplitude of the 
revivals decreases and their width increases as time goes on and the rephasing becomes less complete. Some of these 
features can be understood from a simple analytical approximation introduced in the next section. 

As the coupling strength is increased at a given temperature the behavior of the two-level system becomes increas- 
ingly erratic. Shorter revival times cause successive revivals to overlap and interfere so that the time evolution appears 
irregular. The coupling strength at which the irregularity emerges is closely tied to the temperature: the higher the 
temperature, the smaller the coupling strength needed to produce irregular behavior. Higher temperature also results 
in decreased revival amplitude: larger (N) corresponds to a larger number of frequencies in the sum, which in turn 
causes the rephasing to be less complete. However, a signature of the revivals persists in the form of a return to the 
bare Rabi frequency even at temperatures high enough that the behavior appears random and the revival amplitude 
is essentially washed out. Figure illustrates the lapse into erratic behavior and the persistent revival signature. 

The variation of behavior in the thermal-state adiabatic approximation, from frequency renormalization to distinct 
collapse and revival dynamics to apparent randomness, contrasts with the findings of Knight and Radmore on the 
same type of system within the RWA^ Although they distinguish collapse and revival regions, the behavior within 
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FIG. 7: (Color online) Collapse and first revival for (a) coherent and (b) thermal states in a regime of regular behavior, 
\/loo — 0.1. For large (A'^), the coherent-state model results in a complete collapse, with no residual oscillation amplitude in 
the collapse region; this is not true for the thermal-state model. 



those regions appears erratic, reminiscent of that found above for large coupUng strengths. Clear and well-defined 
revival pulses do not occur in the JCM for any parameter values if the oscillator begins in a thermal state. In fact, the 
basic shape of the time-evolution curve is invariant in the JCM, although the time scale and amplitude may change. 
The distinct revival areas found in the adiabatic approximation for smaller coupling strengths more closely resemble 
those obtained in the coherent-state JCM^^^ 

Finally we turn to the case in which the oscillator begins in a displaced coherent state. Coherent states are considered 
to be the quantum states of the harmonic oscillator which most closely approach the classical limit; it is expected 
that a driven NR is approximated by a coherent state when the external driving dominates thermal fluctuations. The 
initial condition for the coupled system may be written as the pure state 

N=a ^ '' 

where a is the (complex) amplitude of the coherent state, and we will define (N) = |ap. 

The time evolution of the corresponding density matrix Pcoh(O) = |^/'coh(0)) (V'coh(O)l may be calculated in the same 
way as before, with the resulting time-dependent probability 

oo 

i^coh(-, t)^Yl PcohW cos2(f}'Wi/2), (18) 

where PcohiN) — e'~^^'> (N)^ /Nl. Equation 1)18(1 has the same form as Eq. p6(l with the weighting function for the 
thermal distribution replaced by the weighting function for the diagonal elements of the coherent state. 

As before, the sum in Eq. H18|) requires numerical evaluation if no further approximations are made. However, the 
qualitative behavior is more difhcult to classify. The regime (N) w 0.01 behaves much like the same regime in the 
thermal case. Frequency renormalization is the dominant effect, visible for relatively large values of X/loq. 

For values of (N) > 10 and fairly small coupling strengths A, collapses and revivals appear which look similar to 
those found in the JCMiSS The collapses are complete, with virtually no residual amplitude in the collapse region, 
unlike the above results for the thermal state. For small values of (N) the coherent state does not have a wide enough 
spread in frequencies to create a complete collapse, resulting in residual oscillations in the collapse region. Figure 
\7\ compares the coherent-state behavior with the thermal-state behavior. As in the thermal case, a simple analytic 
approximation that will be presented in the next section explains some of these features. 

At larger coupling strengths, however, the coherent-state behavior does not necessarily lapse into irregularity as 
in the thermal case. For large (N), some unexpected results occur as the coupling is increased. The explanation 
for this lies in the non-monotonic dependence of the modified Rabi frequencies fl'^-^^ on both A and N. il'*^^' is 
plotted as a function of X/uiq in Fig. [SJ the values of N illustrated are chosen from the coherent-state distribution 
with (N) = 100 (inset) in order to give a sense of how the spread in frequency corresponds to the distribution in N. 
Figure ini illustrates the resulting behavior of the two-level system in both time and frequency space. The left-hand 
side shows plots of P{—,t) with (N) — 100 for several values of X/loq. Vertical lines in Fig. |S1 correspond to those 
values of X/ojq for which time series are given in Fig. |^ Considering the coherent-state weighting function PcohiN) 
as a function of frequency yields a Fourier-transform- like distribution Pcohi^'^^^), which gives the amplitude of each 
frequency in the sum P(—,t). The weighted frequency distribution corresponding to each time series is plotted on 
the right-hand side of Fig. |51 
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FIG. 8: (Color online) Inset: Normalized coherent-state probability distribution as a function of A'^ with (A*') = 100. Vertical 
lines correspond to the values of whose frequencies are plotted in the main figure. Main figure: Frequency renormalization 
as a function of coupling strength for some representative values of A''. Vertical lines correspond to the coupling strengths used 
in Fig. El 

At zero coupling strength the weighted frequency distribution consists of a delta function located at Q''^^^ /fl — 1. 
As the coupling strength is increased, the center of the distribution shifts toward smaller frequencies and the shape 
spreads out and becomes approximately Gaussian, which results in well-defined collapse and revival regions [Fig.E^a)]. 
When the coupling strength approaches the critical point for the center of the distribution — 0), the function 

Pcoh(f^'^^^) begins to fold back on itself, resulting in a very fast collapse and strangely shaped revivals [Fig.l^b)]. The 
center of the distribution function then shifts back to higher frequencies as the second "hump" of the renormalized 
frequency function fi''^' is traversed (Fig. (S}, and the shape again becomes almost Gaussian. Correspondingly the 
time-dependent probability appears more regular [Fig.E^c)]. At the peak of the second hump of fj'^^^ (A/wq ~ 0.095 
m Fig. El there is very little dispersion in the frequencies corresponding to different values of N, resulting in a nearly 
delta-function distribution for Pcoh (^^''■^'' ) and a very slow collapse [Fig. Efd)]. The frequency distribution function 
then "bounces" back toward low frequencies. Deviations from a Gaussian distribution show up in the time-dependent 
probabilities as altered revival pulse shapes [Fig. I^e)]. The behavior evolves between the examples shown here in a 
continuous manner, although in places it varies quite rapidly with X/uq- 

The previous discussion highlights the complexity of the coherent-state adiabatic model. It should be noted that 
some of the phenomena seen above are blurred out at smaller values of (N). This may be understood by noting that 
two adjacent curves fi'^^) and Jl''^^^^^ diverge more strongly at small values of A'': compare Fig. ^ with N ^ 1, and 
Fig. |S1 with A^ w 100. Thus the weighted frequency distributions Pcoh(^^'''^-') have a wider spread for small values of 
(A'') despite the fact that the coherent-state number distribution is narrower in A^ for small (A^). The wider spread 
in frequency space for smaller (A'^) results in more erratic behavior, without the returns to regularity demonstrated 
in Fig. El 

A similar analysis in frequency space also explains why the thermal-state dynamics appears less complex than the 
coherent-state dynamics. As a function of A^, the maximum of the thermal state distribution is fixed at iV = 
regardless of the value of (N) , while the maximum of the coherent state distribution is given by ( A^ ) . The frequency 
shift in the state A^ = is monotonic as a function of X/ojq and weak compared to the shift for higher values of A^ (Fig. 
21) ■ Thus the maximum of the weighted frequency distribution pth(f^'^^-') for the thermal state shifts less dramatically 
than the maximum of Pcoh(f^'^^-') for the coherent state, and the changes in the shape of the distribution are less 
pronounced. 

Applying the frequency distribution function analysis to the JCM shows how the complicated dependence on A^ and 
A of the frequencies fi'^^' in the adiabatic approximation leads to a much richer variety of behavior than that found 
in the JCM. The Rabi frequencies in the JCM change monotonically with N and A in such a way that the frequency 
distribution function changes in width but never in shape. Changes in (A^) and A in the JCM result in changes in the 
amplitude and time scale of the evolution; however, the shape of the time series is unaltered. The large variation of 
behavior with (A'^) and A is a distinctive feature of the adiabatic approximation, not found in the JCM. 
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FIG. 9: A few samples of the unusual behavior which appears in the coherent-state case with (A'^) = 100 (left-hand side) and 
the associated weighted frequency distributions pcoh(i^''^') (right-hand side). The coupling strengths used are A/ojo = (a) 0.04, 
(b) 0.06, (c) 0.072, (d) 0.096, (e) 0.122. 



IV. THE WEAK COUPLING LIMIT 



Although the adiabatic approximation allows evaluation of the time-dependent behavior of the two-level system, 
the results for thermal or coherent oscillator states are given by infinite sums which must be numerically evaluated. 
This is also true in the JCM, although integral approximation techniques have been used to derive some approximate 
analytic expressions)^'^'* An interesting feature of the adiabatic approximation is that taking the weak coupling limit 
allows closed-form evaluation of the time evolution for both thermal and coherent initial states. 
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FIG. 10: (Color online) Comparison of the first (adiabatic) approximation [Eq. II6I1 I and the second (weak-coupling) approx- 
imation [Eq. 12L)I| 1 for the thermal-state model. Also shown is the envelope function derived from Eq. (I2H by neglecting the 
rapidly-oscillating factor sin(57f). Parameters used are (A'') = 5 and X/luq — 0.1. Although these values are near the limit of 
validity of the approximations used, the agreement up to the first revival is excellent. 

For small values of (A/wo)^ the modified Rabi frequencies given in Eq. 1)12(1 may be approximated 

r!'W«r![l-(iV-}-l/2)(2A/cJo)']. (19) 
This approximation enables the sums in Eqs. H16|l and (|18|) to be carried out Jfi, For the thermal state, Eq. 1)16(1 becomes 



1 lcos{nt)cos[n{2X/ujo)'^t/2] 
th(-, ^) ~ 2 2 1 + A{N){1- 



(1 + 2(A^)) sm{nt) sin[r2(2A/tJo)^V2] 



(20) 



{N))sm'^[n{2X/ujo)H/2] 

A comparison of Eqs. 1(16(1 and 1(20(1 is shown in Fig.^| It is apparent that even for relatively large values of X/ujq the 
approximation to the sum works quite well in the initial collapse region and captures the general qualitative behavior 
of the function. In the limit (N) » 1, Eq. reduces to 



sin Qt 



A{N) siTi[n{2X/LUo)^t/2] 



(21) 



From this formula some of the characteristics discussed in the previous section are immediately evident. Revivals 



occur when fl{2X/uJo)'^t/2 



for m = 1,2, . . ., giving a revival time = 27r/[ri(2A/a;o)^] which decreases with 



increasing X/ljq as expected. The minimum oscillation amplitude in the collapse region scales as 1/{N). However, 
Eq. ((21(1 diverges as i — *■ as well as at the revival times, so it is not useful in predicting the shape of the collapse 
envelope at short times or the nature of the revival envelope function. 

Within the weak-coupling approximation, the coherent-state evolution may be evaluated as well. Equation l(18() 
yields the sum 



Pcoh(-,i)«5 + 



1 -2(A') sin'^[n(2A/wo) t/2] 
2^ 



cos{(iV) sm[n{2X/ujoft] - [1 - {2X/ujQf /2]nt}, 



which in the limit il{2X/u!o)'^t <C 1 reduces to 

Pcoh(-, « i + 1 e-[W"^(2^/-°)'/2l*^ cos{[m + ^){2X/u;of 1]^}- 



(22) 



(23) 



The last form demonstrates the short-time Gaussian collapse envelope which is also found in the coherent-state 
JCMi^iiH?-^'* However, the full expression given in Eq. 1(22(1 is necessary in order to obtain revivals. Notice that the 
revival time is the same as in the thermal case. A comparison of Eqs. ((18(1 and ((22(1 is shown in Fig. ^2 As in the 
thermal state case, the approximation to the sum works well in the initial collapse region. 

This approximation highlights once again the role of the functional form of the modified Rabi frequencies l(12() 
in controlling the time evolution. Rabi frequencies linear in N are obtained in a model similar to the JCM but 
involving three atomic levels and two photons. Since all the frequencies involved arc integer multiples of the bare 
Rabi frequency, the interfering oscillations rephase completely, leading to perfectly periodic, full-amplitude revivals for 
both coherent-state and thermal-state initial conditions. As demonstrated above, the renormalized Rabi frequencies 
in the adiabatic approximation are linear in N to first order. This is not the case in the usual two-level JCM, which 
yields Rabi frequencies which go as the square root of A'^. The closer approach to linearity in N explains why clear 
revivals may be found in the adiabatic model even for a thermal state, while the thermal-state JCM always produces 
erratic behavior. 
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FIG. 11: (Color online) As for Fig. 1101 but for the coherent-state model [Eqs. 11811 and 122ll 1: all parameters are identical, 
envelope shown is the Gaussian collapse function given by Eq. 1231 1 without the rapidly-oscillating cosine factor. 



The 



V. EXPERIMENTAL PROSPECTS 

The primary requirement for experimental exploration of physics in the adiabatic regime is strong coupling at a 
large detuning between the two-level system and the harmonic oscillator. Atom-cavity systems typically have coupling 
strengths A/wq ~ 10"'' - 10"'^ at detunings of A/^o ~ 10"'^ - 10"^, and are well described by the RWA.^'^ The 
adiabatic regime requires numbers several orders of magnitude larger, unlikely to be achieved with atoms. However, 
recent progress in solid-state systems suggests that experimental implementation of the model discussed here may be 
possible fairly soon. A system consisting of a CPB coupled to a superconducting transmission line has achieved the 
"strong-coupling limit" of CQED, in which coherent dynamics occurs faster than the decoherence rates, as confirmed 
by spectroscopio4 and time-domain^ experiments. The coupling strength obtained at zero detuning was A/wq ~ 10~^, 
which is a significant improvement over atomic systems. However, since the coupling is through the electric dipole 
moment it may be difficult to increase the coupling strength by the two orders of magnitude needed to get out of the 
RWA regime, and significant coupling at large detunings is unlikely. 

Coupling strengths much larger than those possible with dipole coupling may be achieved with capacitive or induc- 
tive couplings. Some recent experimental results on a flux qubit inductively coupled to a superconducting quantum 
interference device (SQUID) that acts as both a measuring device and a quantum harmonic oscillator appear very 
promising^ Coherent oscillations in the qubit were observed with Rabi frequencies ranging from very small to as 
large as the qubit splitting frequency. Since the oscillator frequency was about half the qubit splitting frequency, an 
analysis similar to that given below for the CPB-NR system yields values of Q./bjQ < 2. Given the coupling strength 
of X/ujQ w 0.1, the adiabatic regime is already within the reach of this system. 

Although it has not yet been experimentally demonstrated, the system consisting of a CPB capacitively coupled to 
a NR appears to be another potential candidate for achieving the adiabatic regime. The remainder of this section is 
devoted to an analysis of the circumstances under which this would be possible. 

The Hamiltonian is given 



-ffTOTAL = ^^CPB + -f^NR + ^^int 

ffcPB = 4:Ec{ng - \)d, - \Eja^ 
Hint = hX{a^ + a)&^, 



(24) 
(25) 
(26) 
(27) 



where , a are harmonic-oscillator raising and lowering operators which act on the NR; az , d'x are Pauli spin matrices 
operating in the charge basis of the CPB; Ug = {CbVb + CgVg)/2e where Cb and Vb are the CPB biasing capacitance 
and voltage and Cg and Vg are the capacitance and voltage between the NR and the CPB; Ec and Ej are the Coulomb 
and Josephson energies; ujq is the NR oscillator frequency in the absence of coupling; and A = — 4£'(7n^^Aa;zp/fi<i 

NR 
9 

and d is the distance between the NR and the CPB. 

One way of reducing this Hamiltonian to the form of Eq. ^ is to bias the Cooper-pair box to the degeneracy point, 
where Ug — 1/2. However, typical values of Ej /h are on the order of several gigahertz, whereas the highest reported 
nanomechanical resonator frequency is about 1 GHz^ The approximation we have derived in this paper is based on 



where n^^ = CgVg/2e, Axzp = \/ ft/ (2mwo) is the zero-point position uncertainty of the NR with effective mass 
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the assumption that the effective sphtting frequency of the two-level system is much smaller than the frequency of 
the oscillator, f2 ^ wq, which would be difficult to satisfy for degeneracy-point biasing. 

However, an effective Hamiltonian of the correct form which satisfies the condition fl ujq may be found given a 
few reasonable assumptions. Physically, this involves biasing the CPB well away from degeneracy so that AEcing — 
1/2) ^ Ej/2 and applying an oscillating voltage to the CPB bias gate, a procedure which is used in performing 
spectroscopy!^ Assuming that the frequency of the oscillating voltage is equal to the splitting frequency of the box 
and that the amplitude of the oscillating voltage is small, the Hamiltonian may be approximated by Eq. JQ) with 
n = 8gEjEcil/2 - ng)/{[SEcil/2 - Ug)]'^ + Ej}, where g = AEcCbV^^" /2eh is the Rabi frequency induced by the 
oscillating voltage. The Appendix contains a detailed derivation of this approximation. 

Achievable parameters for an experiment of this type align well with the regime in which the adiabatic approximation 
is valid. Starting with a NR frequencjsSiwo/STr = 1 GHz, taking il/wp = 1/10 requires an effective CPB Rabi frequency 
rj/27r = 100 MHz. With CPB parameters Ec/h = 20 GHz, Ej/h = 7 GHz, and Ug = 1/4, the Rabi frequency induced 
by the oscillating voltage should be about 0.5 GHz. Assuming a bias-gate capacitance of Cf, — 10 aP, the required 
amplitude for the oscillating voltage is V^'^ = 0.2 mV, which may be achieved easily. Note that these numbers 
also satisfy the requirements for derivation of the effective Hamiltonian which are given in the Appendix. Since the 
coupling is capacitive, X/loq is limited by how small the distance between the NR and the CPB island can be made 
and by how large a voltage may be applied without damage to the NR. Values on the order of X/ujq w 10^^ - 1 should 
be possible. All of these parameters appear to be well within the reach of present technology. 

Unfortunately the dephasing times for coherent oscillations in a CPB which have been measured so far are quite 
short. Vion et alm^ found a dephasing time of 0.5 ^s at the degeneracy point. As the bias voltage is tuned away 
from the degeneracy point, the dephasing time drops rapidlyji^Si which is attributed to low-frequency charge noise. In 
order to see the effects predicted here, coherence times of several effective Rabi periods, on the order of several tens 
of nanoseconds, would be necessary; this would require significant improvement over current experiments. However, 
charge noise is not believed to be intrinsic to these systems, and advances in materials and fabrication may reduce 
the problem. 

VI. CONCLUSION 

The adiabatic approximation we have discussed in this paper provides a rich and robust framework for exploring spin- 
oscillator physics outside the rotating-wave approximation. Although it is derived under the assumption that the two- 
level splitting frequency is much smaller than the oscillator frequency loq, it works well even when = luq; indeed, it 
provides a more accurate description at large coupling strengths than the RWA. The energy levels obtained from this 
model exhibit nonmonotonic dependence on the oscillator occupation number N and the coupling strength A. This 
leads to complicated time-dependent behavior of the two-level system, which may exhibit frequency renormalization, 
collapse and revival of coherent oscillations, or apparent randomness. Such behavior is quite sensitive to the initial 
state of the oscillator and the coupling strength in some parameter regimes. 

Although the adiabatic approximation is not a new result, solid-state experiments currently underway provide 
motivation for a more thorough exploration of its consequences. The pursuit of quantum computing has catalyzed 
the development of new types of devices which act as artificial atoms. Given the success of atom-cavity experiments 
in demonstrating various characteristics of quantum behavior, it is not surprising that solid-state analogs are being 
pursued. Such systems have the capability to reach regimes, inaccessible to traditional atom-cavity systems, in which 
the RWA is no longer valid. This paper has demonstrated some of the complexity which may be encountered at large 
detuning and strong coupling. 

In particular, we have chosen to focus on a charge-based two-level system coupled to a nanomechanical resonator. 
Observation of the two-level system may offer some insight into the quantum nature of the resonator, just as atoms 
provide a sensitive probe for the nonclassical nature of electromagnetic fields. At fairly high resonator temperatures, 
the shape of the collapse of the coherent oscillations in the CPB may provide some information about the distribution 
of NR states: a thermal state gives a different envelope function than a coherent state. For either distribution, the 
shift from collapse dynamics to frequency renormalization would be a clear indication of near-ground-state cooling of 
the resonator. Finally, the observation of revivals, which are a strictly nonclassical phenomenon, would give evidence 
for the quantum nature of a macroscopic mechanical object. Such experiments appear to be nearly within the reach 
of current technology. 
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APPENDIX: DERIVATION OF THE EFFECTIVE CPB HAMILTONIAN 

The Haniiltonian for the Cooper-pair box in the basis of charge states is^ 

where az,d'x are Pauli spin matrices operating in the charge basis of the CPB; Ug — Cf,Vt/2e where Cb and Vf, are 
the CPB biasing capacitance and voltage; and Ec and Ej are the Coulomb and Josephson energies. Eq. (TOTl may 
also be written in terms of the mixing angle 77 = ta.n~^ {E j / [8Ec {I /2 — ng)]}: 

HcPB = --5A£;(?7)(cos?7(7^ +simjax), (A. 2) 



where lS.E{jf) — ■\/[8£'c'(l/2 — ng)Y + E'^. Alternatively, the Hamiltonian may be written in the diagonal form 

#CPB = -\^E{ri)p\ (A.3) 

where the Pauli operators which operate in the eigenbasis of i/cPB are defined as 

Pz = cos rjcr^ + sin rjCFj, (A.4) 
Px = cos rjcr^ - sin rja^ (A. 5) 

Py = CFy. (A.6) 

An oscillating voltage V^"^ may be applied to the bias gate of the CPB, resulting in an additional term in the 
HamiltonianHi H^^ = Hg cos{u!t) a ^ where g = AEcCtV^ /2eh. In the eigenbasis of i?cPB this becomes 

H.C = i%(e''^* + e-''^*)[cos77p, - sinr;(p+ + /5_)], (A.7) 

where p± = \ {px i ipy) are raising (lowering) operators for the CPB. 

We will assume that uj ~ AE{ri)/h so that the rotating wave approximation may be used to derive a time- 
independent effective Hamiltonian for the CPB with oscillating bias voltage. The first step is to transform into a 
reference frame which rotates about the pz axis at the frequency uj. This may be accomplished by the transformation 
^rot = HU -ihJjUU / dt with U = exp(iwtp^/2). Noting that JJ^^ p± U = exp(=Fit<;t)/5±, the Hamiltonian iJcPB+^^ac 
transforms to 

H,,t = ~\hl\pz + i%(e*"* + e-*"*) C0S77/5, - i%sinr/(p+ + p_ + e-'''^'p+ + e^'^'p^) (A.8) 

where A = AE{ri)/h — w is the detuning between the CPB splitting frequency and the frequency of the oscillating 
bias voltage. As long as we are interested in motion on the timescale of l/g ^ the time-dependent terms in Eq. 
(|A.8p may be neglected. With this approximation, and transforming back to the charge basis, we obtain 

HcPB ~ —^hA{cosri<Tz + siuTyo-j,) — i%sin?7(cos7y(T2. — sinryo-^). (A. 9) 

If we assume the CPB to be biased far from degeneracy such that sin 77 <C cos 77 and take the detuning A = 0, we find 
an effective Hamiltonian for the CPE'^'* 

^CPB = -hf^^^. (A. 10) 

where = g sin r; cos 77. The Hamiltonians for the NR and the interaction between the CPB and NR remain as given 
in Eq. (|24|l . Combining all the terms yields a Hamiltonian for the coupled system of the form of Eq. . 



* Electronic address: |eirish@pas.rochester.edu| 



16 



^ E. T. Jaynes and F. W. Cummings, Proc. IEEE 51, 89 (1963). 

^ J. M. Raimond, M. Bmne, and S. Haroche, Rev. Mod. Phys. 73, 565 (2001). 

3 C. J. Hood, T. W. Lynn, A. C. Doherty, A. S. Parkins, and H. J. Kimble, Science 287, 1447 (2000). 

* A. WallrafF, D. I. Schuster, A. Blais, L. Frunzio, R.-S. Huang, J. Majer, S. Kumar, S. M. Girvin, and R. J. Schoelkopf, 
Nature 431, 162 (2004). 

^ A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, J. Majer, M. H. Devoret, S. M. Girvin, and R. J. Schoelkopf, Phys. Rev. 
Lett. 95, 060501 (2005). 

^ I. Chiorescu, P. Bertet, K. Semba, Y. Nakamura, C. J. P. M. Harmans, and J. E. Mooij, Nature 431, 159 (2004). 

A. D. Armour, M. P. Blencowe, and K. C. Schwab, Physica B 316-317, 406 (2002). 
® A. D. Armour, M. P. Blencowe, and K. C. Schwab, Phys. Rev. Lett. 88, 148301 (2002). 
^ E. K. Irish and K. Schwab, Phys. Rev. B 68, 155311 (2003). 

^° A. T. Sornborger, A. N. Cleland, and M. R. Geller, Phys. Rev. A 70, 052315 (2004). 
" B. W. Shore and P. L. Knight, J. Mod. Optics 40, 1195 (1993). 
R. Graham and M. Hohnerbach, Z. Phys. B 57, 233 (1984). 
R. Graham and M. Hohnerbach, Acta Phys. Austriaca 56, 45 (1984). 
R. Graham and M. Hohnerbach, Phys. Lett. A 101, 61 (1984). 
S. Schweber, Ann. Phys. (N.Y.) 41, 205 (1967). 
M. D. Crisp, Phys. Rev. A 46, 4138 (1992). 

L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambridge University Press, Cambridge, 1995). 
A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. Fisher, A. Garg, and W. Zwerger, Rev. Mod. Phys. 59, 1 (1987). 
D. H. Santamore, A. C. Doherty, and M. C. Cross, Phys. Rev. B 70, 144301 (2004). 
M. D. LaHaye, O. Bun, B. Camarota, and K. C. Schwab, Science 304, 74 (2004). 
F. W. Cummings, Phys. Rev. 140, A1051 (1965). 

J. H. Eberly, N. B. Narozhny, and J. J. Sanchez-Mondragon, Phys. Rev. Lett. 44, 1323 (1980). 
P. L. Knight and P. M. Radmore, Phys. Lett. A 90, 342 (1982). 

N. B. Narozhny, J. J. Sanchez-Mondragon, and J. H. Eberly, Phys. Rev. A 23, 236 (1981). 

M. Abramowitz and I. A. Stegun, eds.. Handbook of Mathematical Functions (Dover Publications, New York, 1970). 
P. L. Knight, Phys. Scr. T12, 51 (1986). 

X. M. H. Huang, C. A. Zorman, M. Mehregany, and M. L. Roukes, Nature (London) 421, 496 (2003). 

K. W. Lehnert, K. Bladh, L. F. Spietz, D. Gunnarsson, D. I. Schuster, P. Delsing, and R. J. Schoelkopf, Phys. Rev. Lett. 
90, 027002 (2003). 

D. Vion, A. Aassime, A. Cottet, P. Joyez, H. Pothier, C. Urbina, D. Esteve, and M. H. Devoret, Science 296, 886 (2002). 
T. Duty, D. Gunnarsson, K. Bladh, and P. Delsing, Phys. Rev. B 69, 140503(R) (2004). 
Y. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. Phys. 73, 357 (2001). 

The same Hamiltonian may be obtained by biasing the CPB at its degeneracy point; however, in this case, typical experi- 
mental parameters result in f2 ^ tjo, while the approximation discussed in this paper works in the opposite limit. See, for 
instance, Ref. |9|. 

It may be interesting to note that 17''"-* may be obtained from Eq. (3.20) of Ref. by setting the spectral density J{uj) = 
5{uj — ujo): that is, by treating the NR as a single-mode bath causing decoherence of the CPB. 

For the sake of simplicity, the effective Hamiltonian is left in the rotating frame. The effect on the expectation value of a^, 
which is all we consider in this paper, is neglible within the assumptions we have already made. 



